Effects of interorbital hopping on orbital fluctuations and metal-insulator transitions: 
Extended linearized dynamical mean-field theory 
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We study the effects of interorbital hopping on orbital fluctuations and Mott-Hubbard metal- 
insulator transition (MIT) in the two-orbital Hubbard model within the extended linearized dynam- 
ical mean-field theory. By mapping the model onto an effective model with different bandwidths 
through the canonical transformation, we find that at half-filling, the increases of the interorbital 
Coulomb interaction U' and the Hund's coupling J drive the MIT, and the critical J c for MIT 
increases with the lift of the inter-orbital hopping integral t a b- Meanwhile at quarter filling and in 
the strong correlation regime, the system without t a \, exhibits MIT with the decreasing of J, and 
favors the orbital liquid ground state. However, the system transits from metal to insulator with the 
increasing of t a (,, accompanied with the rising of the orbital order parameter. These results show 
the important role of the interorbital hopping in the orbital fluctuation and orbital ordering. 

PACS numbers: 71.30.+h, 71.70.Ej, 71.27. +a, 71.20.-b 



I. INTRODUCTION 

In recent years, significant interest has been attracted 
to the problem of orbital ordering and its effect on 
Mott-Hubbard metal-insulator transition (MIT) of some 
transition-metal oxidesi. As a result of the addition of 
orbital degree of freedom, some strongly correlated 3d 
transition-metal oxides exhibit very rich and complicated 
phase diagrams^ii^. In these strongly correlated sys- 
tems with charge, spin and orbital degrees of freedom, 
various orbital, magnetic and charge orders coexist or 
compete with each other. Since many different quantum 
phases are almost degenerate, a little change in orbital 
configuration usually leads to quantum phase transition. 
This has been a hot topic in recent years since the long- 
range orbital ordered ground state is found in colossal 
magnetoresistive parent material LaMnOj2< and the typ- 
ical MIT system The quantum fluctuations are 
important in these compounds, and it comes to be the 
key issue to realize what role the orbital degree of free- 
dom plays in driving the quantum phase transitions. The 
conventional self-consistent mean-field approach under- 
estimates the effect of the orbital quantum fluctuations, 
the dynamical mean- field theory (DMFT) 8 is an useful 
tool to unveil the role of the quantum fluctuations on the 
Mott-Hubbard MIT in these strongly correlated systems. 

The DMFT, which is based on a mapping of lattice 
models onto effective quantum impurity models subject 
to a self-consistency condition, has been very successful 
in describing many aspects of strongly correlated elec- 
tron systems^. Many well-established numerical methods 
have been applied to solve the effective quantum impu- 
rity model, such as the quantum Monte Carlrj2ii2iii, the 
exact diagonalizationi^ii&i^, and the renormalization- 
group method 15,16,17 . Although these methods are suc- 
cessful in studying the MIT in single-band Hubbard 



model, computational limits are met in applying the 
DMFT methods to more complicated multi-band mod- 
els. Recently Bulla and Potthofii^ proposed the lin- 
earized DMFT (LDMFT) to solve the effective impurity 
Anderson model with one impurity site and one bath 
site only, and this ansatz shows minimum computational 
effort in dealing with the MIT of Hubbard model. How- 
ever, the LDMFT is bound to the critical point in pa- 
rameter space. Then Potthoffi^ extended this method 
(ELDMFT) , which is able to access the entire parameter 
space. The ELDMFT reduces to the LDMFT at half- 
filling and at U = U c . The ELDMFT simplifies the pro- 
cedure of DMFT by linearizing the self-consistent equa- 
tions in the low and high energy regions, and can be easily 
applied to study some more complicated models which 
are hard to be treated by the full DMFT. The results 
obtained for the Hubbard model are in good agreement 
with the other numerical techniques^. 

The minimal model to describe the Mott-Hubbard 
MIT including the orbital degree of freedom is the two- 
orbital degenerate Hubbard model. This model has been 
investigated with DMFT by several groupa 20 i 21 i 22 i 23 i 24 i 25 : 
RozenbergSS found that the successive MIT occurs when 
the electron occupation number varies from n=0 to 4; 
Koga et alml showed that two MITs may occur for the 
two subbands at different Coulomb interaction; and Pr- 
uschke et alm^ recently demonstrated that the Hund's 
coupling plays an important role on the MIT. Whereas 
most of these studies neglected the charge transfer be- 
tween the two different orbitals, the inter-orbital charge 
fluctuation was prohibited and the effect of orbital fluc- 
tuation was greatly underestimated. This paper is to 
understand the orbital physics in multi-band Hubbard 
model at finite temperatures, especially for the influences 
of the interorbital hopping and Coulomb interaction on 
the orbital fluctuation, orbital ordering and the MIT in 
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the systems with various electron fillings. 

In what follows, we study the two-orbital Hubbard 
model to explore the role of orbital degrees of freedom 
and the interorbital hopping, Hund's coupling and the 
interorbital Coulomb correlation on MIT both in the half 
and quarter filling cases, respectively. We find the Hund's 
coupling plays different roles in the half- and quarter- 
filling systems; the interorbital Coulomb interaction fur- 
ther splits the two Hubbard bands and drives the system 
from metal to insulator at half filling; the orbital fluctu- 
ation in the presence of interorbital (off-diagonal) orbital 
hopping is so strong that the metallic phase is stable 
at quarter filling; also the interorbital hopping may es- 
tablish the orbital order and drives the system to the 
insulating phase. The paper is organized as follows: in 
Sec. II, we transform the two-orbital Hubbard model with 
interorbital hopping into an effective model with differ- 
ent bandwidth by the canonical transformation, which is 
easily mapped into an impurity model in the framework 
of the DMFT. Then wc describe the extended linearized 
DMFT approach in Sec. III. Our numerical results for the 
systems with half and quarter filling at finite tempera- 
tures are presented in Sec. IV. And Sec.V is devoted for 
the discussion and summary. 



II. THE TWO-ORBITAL HUBBARD MODEL 

The two-orbital Hubbard model consists of the kinetic 
energy and the Coulomb potential parts, 



H = H t + Hu (1) 



with 
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Here I (I') = a or 6, and t aa (tbb) and t a b denote the 
nearest-neighbor intraorbital and interorbital hopping in- 
tegrals, respectively; U describes the on-site intraorbital 
repulsion between electrons and U' the on-site interor- 
bital interaction; and J denotes the Hund's coupling. 
There are two terms with J in Eq. (2). The first one 
corresponds to the spin exchange of the two intraorbital 
electrons, and the second one to the orbital exchange be- 
tween two different orbital electrons. Such a full Hamil- 
tonian keeps the rotation invariance of the Hamiltonian 
(1) in the spin and orbital spaces. Usually the spin ex- 
change parameter is slightly different from the orbital 
one, while for clarity we assume the two parameters are 
the same. For an isolated ion we have the rotation invari- 
ant relation U — U' = 2 J; in some crystals, however, the 



direct Coulomb interactions U and U' are modified by the 
screening effect, and the parameter U' and J should be 
considered as independent parameters^. Here we con- 
sider these two cases, respectively; and for clarity, we 
concentrate on the orbital fluctuation, and fix the spin 
configuration as the paramagnetic phase, that is, we look 
for dynamic mean-field solutions that do not break the 
spin rotational symmetry in this paper. 

For such multi-orbital systems with off-diagonal orbital 
hopping, it is hard to map the Hamiltonian (1) onto 
the single impurity effective model. We introduce two 
fermions a' and fy via a canonical transformation, and 
the kinetic terms become of orbital diagonal in this quasi- 
particle representation. The canonical transformation is 
expressed as: 



Ciotr = ua ir7 + vf3 itT 
Ciba = ~va ia + uf3 ia , 



(3) 



with 




t 2 ab + (*oa - tbb) 2 



(4) 



Here a 1 and j3' are the new fermion operators, respec- 
tively. In this situation the twofold degenerate orbital 
degree of freedom is converted into two different kinds of 
fermions, and the orbital correlation is transformed into 
the particle-particle correlation. By Eq.(3), the kinetic 
energy H t in Eq.(l) is expressed as, 



j & 



<i,j>,<7 <i,j>,<7 



with t a = t aa u 2 + tbbV 2 — t a buv and tp = t a 



t b bU 2 



t a buv. In addition, we transform the potential energy 
Hjj via the canonical transformation Eq.(3), and ob- 
tain an effective two-orbital Hamiltonian with different 
bandwidths so long as t a b ^ 0. Through the canonical 
transformation the two-orbital fermions with interorbital 
hopping are transferred into two new fermions without 
interorbital hopping, which can be treated more easily 
in the DMFT scheme as to the single Hubbard model. 
On the other hand, the on-site Coulomb energy and the 
Hund's coupling in Hjj inevitably become very compli- 
cated, which is not hard to treat in the DMFT framework 
in the present fermion representation. 
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III. THE EXTENDED LINEARIZED DMFT 

In the DMFT, the above effective lattice model can be 
mapped into an impurity model with two orbitals, 



H, 



imp — /] eiiaa>fa rT a h(T ~ M / ^{ a og Q! Og + PoaPov} 

+ J2^°Ma^0a + h.c) + V^(a+j0a + k.c)} 



(6) 



where ao a (or /3o<t) annihilates a fermion with spin a at 
the impurity site i = 0. The impurity couples with the 
bath which is described by operators (a^, a; 7<J ) with 
energy e/ 7CT via hybridization V^ a (7 = a or 0). For this 
model, the free (U = U' = J = 0) impurity Green's 
function G°(7,7') = {{joa, lot)) (loa = a 0a or 

/3o<r) is a 2 x 2 matrix, and it can be obtained via, 



(G^M)- 1 







m-E 



(vj'J 2 

la u)-t Wa 



(7) 



The Green's function G a (uj) of the impurity model in 
Eq.(6) is obtained by the exact diagonalization method 8 , 
and thus the self-energy matrix is directly extracted by 
Dyson equation, 



s (T M = (G°H)- 1 -(G CT H)- 1 . 

Therefore, the on-site Green's function reads, 

G.M = r de P0ie) y(lY 
J -00 w + /i-e-E CT (cj) 



(8) 



(9) 



Here we choose the semicircle density of states (DOS) 
Po(e)a,p = \J^t 2 a .p — e2 /(27r^ i /3) with respect to the two 
branch fermions. The self-consistent equations (8) and 
(9) are satisfied for proper values of the parameters ei ia 
and V^ a of the impurity model by iterative calculation. 

The exact diagonalization method has been combined 
with DMFT, in which finite orbitals are chosen to de- 
scribe the bath&. This method is exactly only as the 
number of the bath orbitals Nb=oo. However, in the 
study of single band Hubbard model, it has been found 
that when Nf, > 5, the corresponding results are good 
enough. While for more complicated models, such as the 
multi-orbital degenerate models or large cluster "impu- 
rity" model, the corresponding orbital parameter space is 
much larger than that of the single-band Hubbard model, 
and the self-consistent iteration will consume huge com- 
putation time. In order to obtain precise results within 
the ability of implementations of DMFT, the ELDMFT 
method is constructed by considering an effective im- 
purity model with one impurity site and one bath site 
onljiiS. This method simplifies the procedure of DMFT 



by linearizing the self-consistent equations in the low and 
high energy regions, and two self-consistent conditions 
are introduced to fix the two bath parameters e and V. 
For example in the single-band Hubbard model, Potthoff 
compared the high-frequency expansion of the on-site lat- 
tice Green's function with the result of the exact high- 
frequency expansion, and found that the electron fillings 
given by the Green's functions and by the impurity model 



is equa 
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(10) 



Eq. (10) is the first self-consistent equation for the EL- 
DMFT. To satisfy the DMFT self-consistent equation 
(Eq. (9)) in the low frequencies and match the features 
of the central quasiparticle peak, Potthoff introduced the 
second self-consistent equationi^, 



V 2 = ZM. 



(«) 



(11) 



with M 2 ° = J dee 2 po(e). For a metal, Z denotes the 
quasi-particle weight, which represents the single parti- 
cle excitation near the Fermi surface of the metal. As Z 
approaches zero, the single particle excitation vanishes, 
indicating the occurrence of MIT. The results for the 
single-band Hubbard model are in good agreement with 
other numerical techniques 19 . 

To study the two-orbital Hubbard model, we introduce 
four bath parameters (e a , ep, V a , and V 13 ) to describe 
the paramagnetic phase of the bath, and thus four self- 
consistent conditions are introduced to fix them accord- 
ingly, 



i 7 = 

imp 



Z^M, 



(0) 



(12) 



with — J dee 2 po(e) 7 , and 7 = a or (3. These equa- 
tions renormalize the high- and low-frequency regimes 
of the original self-consistency conditions in an integral, 
qualitative form, and are thus well motivated. 



IV. THEORETICAL RESULTS 

In studying the orbital fluctuations and MIT of the 
two-orbital Hubbard model at finite temperature by the 
ELDMFT method, we choose the original diagonal hop- 
ping as the unit, t aa = ^fcb=l, and mainly consider two 
different conditions: one is the spin symmetric case with 
2, J = U — U', and the other one is the case with the 
parameters J and U' being independent. Without spe- 
cial note being made, the temperature in this paper is 
fcgT = 1//? = 1/16. We restrict the present study to the 
paramagnetic phase, the other possible magnetic ordered 
phases shall be encountered latter on. In the following 
studies, we focus on the systems with half filling (ri = 2) 
and quarter filling (n = 1), respectively. 
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FIG. 1: Evolution of the DOS with the interorbital Coulomb 
interaction in the half-filling systems. Theoretical parameters 
are {7=4, J = and t a b = 0. 



A. Half filling cases 

The two-orbital Hubbard model reduces to two inde- 
pendent single-band Hubbard models at t a b = U' = J = 
0. For the symmetric case with half filling n = 2, the MIT 
is observed at the critical Coulomb interaction U c w 6.0, 
which is consistent with the Potthoff's result. We now 
consider the systems with t a b = J = and U = 4, and 
study the effect of the interorbital Coulomb interaction 
on the MIT and orbital fluctuations. 

The DOS for the systems with U'=0.0, 2.0, 4.0, 4.7 
and 5.0 are shown in Fig. 2] At U'=0 we find that the 
quasiparticle spectra exhibit three distinct features: the 
broad upper and lower Hubbard bands and the dominant 
quasiparticle band near the Fermi energy, indicating that 
the system is a good metal. This result resembles to that 
of the single-band Hubbard model by the full DMFT&. 
For the systems with finite U', the spectral characters 
in the regime of V < U' c is quite different from that in 
V > U' c , here U' c denotes the MIT critical point. In the 
regime of U' < U' c , our results clearly demonstrate that 
the interorbital correlations remove the orbital degener- 
acy of the Hubbard bands and further split the upper and 
lower Hubbard bands into four Hubbard subbands. Ap- 
proaching the critical point U' c , the DOS near the Fermi 
surface decreases very quickly with the increase of U', as 
shown in Fig. Large U' suppresses the quasiparticle 
excitation near E p , and as a result it leads to the MIT at 
the critical point of U' c = 4.75. What we obtained is in 
agreement with the previous results by RozenbergSS and 
Koga et alM>, confirming the validity and the reliability 
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FIG. 2: Quasiparticle weights Z (dotted line), Z a (filled sym- 
bols) and Zp (open symbols) as the function of U' at half 
filling with different interorbital hoping t a t- The other pa- 
rameters are: U=4 and J=0. 



of the present ELDMFT method. 

To learn the effects of the orbital fluctuations more bet- 
ter, we also study the influence of thermal fluctuation on 
MIT at low temperature /3 = 32 and high temperature 
/3 = 8 by keeping the other parameters fixed. Consid- 
erable influence of the thermal fluctuations are found, 
and the corresponding upper critical points of MIT are 
[7,1=5.0 for (3 = 32 and [7^=4.45 for /3 = 8, respectively. 
We do not consider a lower critical point of MIT with 
U' > U, presumptively unphysical. Obviously, the in- 
terorbital interaction drives the MIT at high temperature 
more easier than at low temperature, since the combina- 
tion of the quantum and the thermal fluctuations drives 
the paramagnetic metal transition to the paramagnetic 
insulator more easily. 

Next we study the influence of interorbital hopping on 
the orbital fluctuations and the MIT in the system with 
half filling. The dependences of the quasiparticle weights 
on the interorbital interaction in the cases with t a b = 0.5 
and tab = 1.0 are shown in Fig. [5] In Fig. [2J the two 
dashed lines with open and filled circles represent the 
quasiparticle weight Z of the effective orbitals a and f3 
as t a b = 0.5 respectively, while the two solid lines with 
open and filled triangles arc the corresponding results for 
tab = 1-0. For the convenience of comparison, the result 
of the degenerate bands as t a b — (dotted line) is also 
shown. Firstly, we find that the difference between the 
two effective bands a and (3 is significant as the interor- 
bital hopping approaches the intraorbital hopping. In 
comparison, the quasiparticle weights of the a- and fi- 
orbitals are identical for the system without off-diagonal 
hopping. In the presence of the interorbital hopping, the 
quasiparticle weight of the /3-band is larger than that of 
the system with t a b = in the metallic phase, while the 
a-band is smaller. As shown in Fig. [21 the contribution 
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FIG. 3: (a) Dependence of quasiparticle weights Z (dotted 
line), Z a (filled symbols) and Zp (open symbols) as the func- 
tion of J at t ao = 1-0, and (b) the critical value J c as a func- 
tion of interorbital hopping integral t a b in the systems with 
2J = U-U'. 



of the a-band to the metallic properties at t a t — 1.0 is 
small, the properties of the system is mainly controlled 
by the /3-band, implying that the orbital excitation and 
orbital fluctuation are large and dominant in the metal- 
lic phase. Secondly, the MIT critical point of the sys- 
tem is U' c = 4.5 for t a b = 1.0, obviously smaller than 
U' c = 4.7 for t a b = 0.5. The interorbital hopping broad- 
ens the bandwidth of /3-band and narrows that of a-band, 
and more electrons transfer to the /3-band with the in- 
crease of t a b. As we have obtained above, the interor- 
bital interaction suppresses the quasiparticle excitation 
near Fermi surface in the cases without off-diagonal hop- 
ping. Similarly, the interorbital hopping also depresses 
the quasiparticle weight of the a-band by narrowing its 
bandwidth. As a result, the MIT is more easier to happen 
when t a i, becomes large in the paramagnetic half-filling 
system. 

Apart from the effects of interorbital hopping and 
Coulomb correlation on MIT, our calculations also show 
that the Hund's coupling J plays a very important role 
in MIT. We find that for t a b—0, the behaviors of Z for 
the J ^0 case are quilt different from the J = case. 
There are two MIT critical points with J c = 0.55 and 
J' c = —0.25 (unphysical) for the spin rotation symmetry 
case 2 J = U—U', and the metallic phase can exist only in 
a narrow region J' c < J < J c , which is in agreement with 
Pruschke and Bulla's result. The large Hund's coupling 
J suppresses spin flip excitation, and as J increases the 
exchange splitting pushes the a-bands to higher with re- 
spect to the /3-bands, resulting in the system to transition 
to the insulating phase. Therefore the role of the Hund's 
coupling is similar to the on-site Coulomb Interaction U , 
driving the MIT by preventing the electrons from double 
occupation. Meanwhile, instead of excluding the double 
occupation, the interorbital Coulomb correlation drives 
the MIT by removing the orbital degeneracy of the two 
Hubbard bands. 



Further more we study the interplay between the in- 
terorbital hopping and the Hund's coupling in driving 
the MIT. We show the dependence of the quasiparti- 
cle weight on the Hund's coupling in the system with 
2J = U - U' and t ab = 1.0 in Fig. OJa). For the conve- 
nience of comparison, the result of the degenerate bands 
as tab — (dotted line) is also shown. When J is small, 
the system is metallic, and the quasiparticle weight of 
the effective /3-band is considerable larger than that of 
tab = 0. Since the hopping channels of the former is more 
than the latter, the charge transfer between two different 
orbitals leads to more considerable orbital quantum fluc- 
tuation than that in t ao = 0. That means the presence of 
t a b broads the conduction bandwidth and increases the 
DOS of the /3-band, and thus enhances the quasiparticle 
weight Z@, as shown in Fig. 2. When J becomes large, a 
lot of spin flip excitation is suppressed, the system under- 
goes a MIT, correspondingly defines a critical value J c . 
The MIT critical point J c as a function of the interorbital 
hopping is shown in Fig. |3{b). Obviously the increase of 
the interorbital hopping enhance the orbital fluctuations, 
and contributes more excitations near the Fermi surface. 
As a result, more large J c are needed to drive the MIT, 
and thus the MIT critical value J c increases with the lift 
of the interorbital hopping. 



B. Quarter filling cases 

The two-orbital Hubbard model with quarter filling is 
appropriate to describe some important strongly corre- 
lated materials with twofold degenerate orbital degree 
of freedom, such as the hole-type compound KC11F3 and 
electron-type compound LaNiO^. To our knowledge, the 
strongly correlated systems with quarter filling have not 
been given enough attention by the conventional DMFT 
yet. In this subsection we study the two-orbital Hubbard 
model with quarter filling, we take large intraorbital in- 
teractions (U — 8 and 12), and assume the relationship 
U = U' + 2J. 

First we consider the case with t a b = 0. The evolution 
of the DOS with the Hund's coupling J {J — 0.0, 0.5, 1.0 
and 2.0) for U—% and 12 are shown in Fig.0] The Hund's 
coupling strongly modifies the orbital subbands. At J=0, 
the presence of the interorbital and intraorbital Coulomb 
interactions U' and U removes nearly a half of the quasi- 
particle weight in the lower Hubbard band to the upper 
Hubbard band, and the lowest subband is thus almost 
filled up and the residue quasiparticle weight near the 
Fermi surface is very small. Therefore the system with 
J = at quarter filling is a bad metal or an insulator. 
In the system with U—8, the DOS at the Fermi level 
increases gradually as J increases; as J > 0.5 the system 
transits from bad metallic or insulating phase to metallic 
phase. Meanwhile, the increasing of Hund's coupling can 
further split the energy bands into many subbands as 
shown in Fig.^ It can be seen clearly that in the quarter- 
filling case, the influence of the Hund's coupling J on 
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FIG. 4: Evolution of the DOS with the Hund's coupling J 
the quarter-filling systems with U — 8 and U — 12 at t a b=( 
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MIT is different to the half-filling case: in the former tl 
increase of the Hund's coupling indicates the reduction 
the interorbital repulsion U' and the enhancement of the 
ferromagnetic fluctuation, therefore the metallic phase is 
in favorable; however in the latter, the increase of J leads 
the system to transit from metal to insulator since the 
spin flip excitation is suppressed and the quasiparticle 
weight near Fermi surface greatly declines. On the other 
hand, for larger intraorbital Coulomb interaction U = 12, 
the system is always an insulator with the increasing of 
J up to 2. As J approaches to 2, the edge of the upper 
Hubbard bands is gradually shifted to the Fermi energy, 
indicating the MIT will occur at J c w 2. 

From Fig. 4 one can see the striking difference of the 
quasiparticle weights between the systems with U = 8 
and with U—12. At [7 = 8, the quasiparticle weight Z 
monotonously increases with the lift of </, implying that 
the system is a good metal in the absence of the interor- 
bital hopping. However, for the system with U = 12, 
the quasiparticle weight Z is very small, and the total 
weight is 0.031 as the Hund's coupling J increases to 2. 
The double occupancy D, as one expects, decreases with 
the increasing of the intraorbital repulsion. At U = 8, the 
maximum magnitude of the double occupancy is about 
D m = 0.0012. While at U = 12 the double occupancy 
D m , about 0.0001, is negligible. Such strong single occu- 
pation and large Coulomb interaction in the System with 
quarter filling suggests that the system with [7=12 is in 
the antiferromagnetic insulating phase, while the system 
with U =8 is a paramagnetic or ferromagnetic metal. We 
will discuss the cases with spin symmetry broken in a 



future paper. 

In the two-orbital Hubbard model at quarter filling, it 
is very interesting how the interorbital hopping affects 
the MIT. Here wc mainly study the MIT in the systems 
with [7=8 and J=1.0. The evolution of the DOS of the 
a- and /3-bands with the interorbital hopping is shown in 
Fig- At small t a b both a- and /?-bands cross the Fermi 
surface, indicating that the system is metallic. Signif- 
icantly different from the half-filling cases, the further 
increase of the interorbital hopping greatly reduces the 
quasiparticle states near the Fermi surface for both the a- 
and /3-bands. One finds that in the a-band most states 
move to the high energy regime above the Fermi level, 
while in the /3-band, most states move to the low en- 
ergy regime below the Fermi level, thus the asymmetry 
between the a- and the /3-band appears. And we find 
that as tab > 1-0, the DOS of the a- and /3-bands almost 
vanish at the Fermi surface, implying an energy gap will 
open for the two bands and MIT will occur. Since the 
electron kinetic energy increases with t ao , it is out of 
the expectation that the rise of t a b leads to MIT. As we 
show in the following, this arises from the formation of 
the orbital ordering in quarter-filling systems with large 
interorbital hopping integrals. 

In the presence of the strong Coulomb correlation, the 
system with symmetric orbitals tends to break the orbital 
symmetry to lower the ground state energy and form 
long-range orbital order. The orbital order parameter, 
< T >, is usually defined as the average of the pseu- 
dospin operator T = ^ ab Cj a a a bCib, with a being the 
Pauli matrix, representing the orbital polarization of the 
electrons occupation in the two orbitals. In the quarter- 



7 



.10 




FIG. 6: Dependence of orbital order parameter T x on the 
off-diagonal orbital hopping integral t a b in the quarter-filling 
systems with U = 12. 



filling system with U=12, we calculate the three compo- 
nents of T, T x , T v and T z , and the results are shown in 
Fig. From the preceding study, we have found that in 
the absence of the interorbital hopping the system with 
large interorbital interaction is in a para-orbital phase 27 . 
Switching on the off-diagonal or interorbital hopping, we 
find that the symmetry of the orbital space is broken, 
and the x-component of orbital order parameter T x is 
finite, while the z-component T z is always zero. This 
result is also confirmed by the self-consistent mean-field 
method for the effective spin-orbital superexchange in- 
teraction in the limit of large U2£. In addition, large in- 
terorbital hopping results in strong orbital fluctuations, 
and large Hund's coupling favors the ferro-orbital occu- 
pation, the competition of these factors favors the orbital 
polarization in the x-direction, leads to the orbital order 
parameter T lying in the x-direction. and T x increases 
with the increasing of the interorbital hopping integrals 
and the Hund's coupling J, as seen in Fig. |S| 

In the absence of the interorbital hopping t a b, the in- 
terorbital Coulomb interaction U' and the Hund's cou- 
pling J adjust the relative positions of the four subbands; 
though the electrons occupying one orbital configuration 
is favorable in energy, no electron can transfer between 
these orbitals, and the orbital occupations of the elec- 
trons in the two orbital are the same, giving rising to 
orbital disordered phase. In the presence of t a b, the 
electrons can transfer between orbitals and occupy the 
most favorable orbital configuration in energy, and there- 
fore the orbital ordering can establish, as we find in the 
present results. It is worthy of noting that we only con- 
sider a possible ferro-orbital phase, while more compli- 
cated orderings are out of reach from the present single- 
impurity DMFT study. Further study will be addressed 
in the future paper. 



V. CONCLUSIONS 

From the proceeding study it is clear that the interor- 
bital hopping plays an important role in the MIT of the 
two-orbital degenerate Hubbard models. The presence 
of the interorbital hopping not only enhances the orbital 
fluctuations and quasi-particle excitation near Fermi sur- 
face, but also leads to some unusual results. At half 
filling, since more and more electrons transfer from one 
orbital to the another, the critical value of the interor- 
bital Coulomb interaction U' c for MIT decreases with the 
increase of the interorbital hopping; while the critical 
value of the Hund's coupling J c for MIT lifts with the in- 
crease of tab- In the systems with 2 J = U — U' , metallic 
phase can exist in a very narrow region with the vari- 
ation of J, since both intra- and interorbital Coulomb 
interactions, U and U' , split the Hubbard subbands and 
drive the MIT. At quarter filling, the increase of J may 
lead to insulator-metal transition, and this role in MIT 
is completely contrary to that in the half-filling case; and 
the increase of the interorbital hopping also may drive 
the MIT due to strong orbital fluctuation, and leads to 
weak ferro-orbital ordering in the system with large on- 
site Coulomb interaction. 

Meanwhile, the interorbital Coulomb correlation also 
drives the MIT by expanding the separation of two or- 
bital Hubbard subbands both in half- and quarter-filling 
systems. In the present paper we mainly concentrate our 
study on the paramagnetic and para/ferro-orbital phase, 
little attention is paid to complicated spin and orbital or- 
derings. At present, it is not possible to compare the the- 
oretical results with the experimental results in candidate 
compounds KCuF 3 and LaNi03 , since more complicated 
spin and orbital orders are involved in these two realistic 
compounds, thus more than two impurities must be con- 
sidered in the study, beyond the single-impurity DMFT 
theory. 

Recently it receives great interest whether the wide and 
the narrow Hubbard subbands exhibit MIT separately, 
i. e. the orbital-selective Mott transition (OSMT) 22 i 28 in 
two-orbital Hubbard model. Obviously our results do not 
show that there exist such Mott transition in the pres- 
ence of interorbital hoping t a b ^ 0. While in the absence 
of t a b, it is not difficult to expect the existence of OSMT 
since the charge transfer between orbitals is forbidden, 
and the role of the a-band is to exert an effective poten- 
tial on the 6-band and to change the chemical potential 
of the subband, therefore the MIT of the a-band and 
the 6-band occur separately. On the contrary for t a b 
there exists charge transfer between the two orbitals, and 
the MIT will occur at the same critical points for both 
bands. While for the system involving two subbands of 
different bandwidths and omitting the interorbital hop- 
ping, as we find in the new quasi-particle representation 
in Fig. 3a, the OSMT seems happen as the Hund's cou- 
pling is large enough. Further studies are deserved to 
investigate the OSMT in the whole parameters space to 
determine the phase diagram of the two-band Hubbard 
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model with nondegenerate subbands. 

In summary, utilizing the extended linearized DMFT, 
we find in the two-orbital Hubbard models with the in- 
terorbital hopping and half filling, the increase of Hund's 
coupling drives the MIT; on the other hand, the quarter- 
filling two-orbital systems remain metallic due to large 
orbital fluctuations. The systems with only intraorbital 
hopping favors metallic and orbital liquid phase; as a 
contrast, there exists the long-range orbital ordering in 
the quarter-filling two-orbital systems with interorbital 



hopping and large U. 
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